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ABSTRACT 

Stationary solutions of spherically symmetric accretion processes have been subjected to a 
time-dependent radial perturbation, whose equation includes nonlinearity to any arbitrary or- 
der. Regardless of the order of nonlinearity, the equation of the perturbation bears a form that 
is remarkably similar to the metric equation of an analogue acoustic black hole. Casting the 
perturbation as a standing wave and maintaining nonlinearity in it up to the second order, 
bring out the time-dependence of the perturbation in the form of a Lienard system. A dy- 
namical systems analysis of this Lienard system reveals a saddle point in real time, with the 
implication that instabilities will develop in the accreting system when the perturbation is ex- 
tended into the nonlinear regime. The instability of initial subsonic states may also adversely 
affect the temporal evolution of the flow towards a final and stable transonic state. 
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1 INTRODUCTION 

In the context of astrophysical fluid flows, the classical model of spherically symmetric accretion, proposed bv iBondj 1 1952 ) sixty year 
ago, is in essence a mathematical problem of conservative and compressible hydrodynamics. This model has acquired the status of a 
paradigm in studies on accretion, and apart from the fact that it is amenable to exact mathematical analyses in many cases, the spheri- 
cally symmetric model faithfully captures much of the physics of many astrophysical flows. So, nothwithstanding its apparent simplic- 
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ity, the spherically symmetric flow has been a subject of enduring interest to the researcher in astrophysical fluid dynamics, with mul- 
tiple physical and mathematical variations on the original theme ( Park el 1 19581: ISalpetej 1 1964 |Parkej|l966l: Uxford & Newmaij|l967l : 

ian|[l978 



Cowie et al.1 



Holzer & Axford 197G; Balazs 1972; Michel 1972; Meszaros 1975; Blumenthal & Mathews 1976; Meszaros & Silk 1977; Begelm; 

Garlickll 19791: 



19781; Istellingwerf & Bufdl 19781; 



1984; BonazzolaetalJll98 



7l 1 1992l:l 



Theuns&Davidfl 19921- 



); McC ra^ll979l;lBriiikmanrJll98d:lM^^ et al.lll980; 

; Kazhdan & MurzimJI 1994 iRufferJI 1994 lMarkovidll995l: iTsuribe et al 



Titarchuk et al. 1996; Zampieri et al. 1996; Titarchuk et al. 1997; Kovalenko & Eremin 1998: 



Vitellc 



1W4; Markovic I 1 
Da j| 19991: iMalepfT; 



1999; Toropin et al 



1995 



1999 



Dasll200d:lDas & Sarkaill200ll:lRav & Bhattachariej|2002l:lFoglizzcjl2002l:lRavll2003l:lBabichev et alJl2004lDasll2004lRay & Bhattachariee 
20051: lGaitJ200d:lMandal et al.ll2007l:lRovll2007l : lRov & RavlkoollNaskar et alJbool Isilich et alJl2008l :lMach & Maledl2008l : lRovll201ll : 



Park & Ricottill201 lklWong et al.11201 lh . 

Accretion processes involve the flow dynamics of astrophysical matte r under the extern al gravitational influence of an astrophysical 
object, like an ordinary star or a white dwarf or a neutron star or a black hole jFrank et al J2002ri . Accretion flows are distinctly different from 
the self-gravity driven collapse of a fluid system, such as a star. The accreting astrophysical matter could be the interstellar matter, as modelled 
by its spherically symmetric infall onto an isolated accretor, or stellar matter, as seen in a binary sy stem, where the tidal deformation of a 
star causes matter to flow out of it into the potential well of a compact companion {Frank et al]2002f) . In all of these cases, the mathematical 
description of the fluid syst em involves a mo mentum balance equation (with gravity as an external force), the continuity equation and a 
polytropic equation of state jFrank et al ]2002o . 

Fluid flows, conservative or dissipative, fall under the general class of nonlinear dynamics. Set in full detail, the condition for mo- 
mentum craisen^tionjna fluid is a balanc e of dynamic effects, nonlinear effects and the effects of the pressure inherent in a continuum 
system 1 Landau & Lif shitz 1987 ). Prior to Bondil 1952h . some studies of astrophysical flows had considered only the interplay between 
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dynamics and nonlinearity, neglecting the effects of pressure with the argument that any heat generated would be radiated away rapidly, so 
that the temperature of the infalling gas (and relate d to it, the pr essure as well) would remain negligibly low. The other extreme of ignoring 
dynamics and introducing pressure was taken up by Bondi ( 1952), in what became a stationary mathematical problem. Nonlinearity, however, 
was an abiding presence in either case. 

While solving the stationary, spherically symmetric, compressible fluid flow was not difficult mathematically, interpreting the behaviour 
of the solutions from a physical perspective was. From the plethora of mathematical solutions in the stationary problem, the ones of physical 
relevance, in view of modelling astrophysical inflows, were identified to be locally subsonic very far away from the accretor. Within the class 
of inflows obeying this outer boundary condition, there is an infinitude of globally subsonic solutions, along which a fluid element may reach 
the accretor with a low subsonic velocity. However, for the same outer boundary condition, a single critical solution stands out in a class by 
itself, capable of allowing matter to reach the accretor with a high supersonic velocity, and crossing the sonic horizon along the way. This is 
the unique transonic solution — the classical Bondi ( 1952) accretion solution. 



The exact fashion in which accreting matter reaches the accretor is relat ed to the inner b oundary c ondition of the inflow proble m. In the 
case of the accretor being a black hole, the infall process must be transonic I INovikov & Thornej|l973l ; IShapiro & Teu kolsky 1983). This is 
because a black hole has an event horizon instead of a physical surface, and thus precludes all possibility of a pressure build-up at small radii, 
that could otherwise have dominated over the free-fall conditions close to the accretor. The situation, however, is not so clearly understood if 
the accretor has a hard surface like a neutron star or a white dwarf. For such an accretor, it is supposed that the accumulated matter would build 
up pressure near the surface and cause the supersonic flow to be shocked down to subsonic levels, although for a neutron st ar in particular, 



all ac creted matter is expected to be efficiently "vacuum cleaned" away, making it easier for the flow to remain supersonic dPetterson et al 



1980). Evidently then, questions of setting the inner boundary condition an d determining an inflow trajectory in relation to it, are not trivial 
ones to confront. Nevertheless, working with the stationary problem itself. iBondil dl952l) had the insight that the transonic solution would 
be the one selected by a fluid element to reach the accretor from a distant outer boundary. The governing principles behind thi s choice were 
connected to the maximisation of the mass accretion rate, and the minimisation of the total energy configur ation of the flo w ((Bondi 1952; 
iGarlickll 19791) , although a definitive conclusion regarding the realisability of the transonic solution was left bv lBondil d 19521) to its stability. 

The problem with the transonic solution in the stationary regime is that its realisa bility is notoriously vulnera ble to even an infinites- 
imal deviation from the precisely needed boundary condition to generate the solution jRav & Bhattachariee 2002 ). This difficulty may be 
overcome by looking at the possibility of a temporal evolution of the accreting system towards the transonic state 1 Ray & Bhattacharied 



20021; IRov &"Ray]|2007h . However, the nonlinear equations governing the temporal evolution of the flow do not lend themselves to ready 



mathematical analyses. Indeed, in the matter of inc orporating both the dynamic and the pressure effects in the equations, the mathematical 
problem was very aptly described by Bondi ( 1952) as "insuperable". So in the absence of any analytical formulation of the dynamics of 
the flow sol utions, much of all time-dependent studies in spherically symmetric accretion i s perturbative in character, based on linear s t abil- 
ity analysis Jstellingwerf & Bufjl 19781 : lGarliclj[l979l ; betterson et aill980l ; iRuffertll 19941 ; iKovalenko & Erermnlfl^j; iFoglizzolbooj iRav 
2003; Gaite 2006; Roy & Ray 2007), although in this respect, some non-perturbative studies have also been reported JRav & Bhattachariee 



2002 



Roy & Ravll2007h . The range of the perturbative studies covers numerical and analytical methods, using both radial and non-radial 



perturbations, leading to varied conclusions about the behaviour of the perturbations in the spatial and temporal domains. The commonly 
accepted view to have emerged from all the linear stability analyses is that perturbations on the flow do not produce any linear mode with an 
amplitude that gets amp lified in time 1 Gaitdboolil) . and that the perturbative method does not indicate the primacy of any particular class of 
solutions dGarlickll 19791) . This is as far as one could say, working in the linear regime. However, the general experience associated with any 
nonlinear system (and accreting systems are very much nonlinear) is that the understanding gained through linearised conditions can scarcely 
be imposed on circumstances dominated by nonlinearity. The work presented in this paper makes an attempt to br idge this gap. 

In this work, a time-dependent, radial perturbation scheme implemented originally by Petterson et al. 1 1980h has been adopted and all 
orders of nonlinearity have been retained in the resulting equation of perturbation. A most striking feature of the equation of the perturbation 
is that even on accommodating nonlinearity in full order, it conforms to the structure of the metric equation of a scalar field in Lorentzian 
geometry. This fluid analogue (an "acoustic black hole" ), emulating many features of a general relativistic black hole, is a matter of continuin 
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interest in fluid mechanics from diverse points of view (Moncrief 1980; Visser 1998; Schutzhold & Unruh 2002; Barcelo et al. 2005; Volovi 
20051: Isingha et alfeooj iRav & Bhattacharieell2007allbl;lRov & Ray||2007l;lNaskar et alJl2007l ; lDas et alj2007l ; lMach & MalecM2008l) . 



The equation of the perturbation is then applied to study the stability of globally subsonic stationary solutions. Regarding the non- 
perturbative evolution of the accreting system, it is feasible to suggest that the initial condition of the evolution is a globally subsonic state, 
with gravity subsequently driving the system to a transonic state, sweeping through an infinitude of intermediate subsonic states. So, to 
ensure an unhindered temporal convergence to a stable transonic trajectory, the stability of the subsonic states is essential. To investigate this 
aspect at a relatively simple level, all orders of nonlinearity beyond the second order have been truncated in the equation of the perturbation. 
Following this, the spat ial dependence of th e perturbation has been integrated out with the help of well-defined boundary conditions on 
globally subsonic flows JPetterson et al .Il980l) . After this, only t he time-dependent part of the perturb ation is extracted, and, very intriguingly, 
it acquires the mathematical appearance of a Lienard system dStrogatz 1994; Jo rdan & Smithlll999l) . Application of the common analytical 
tools of dynamical systems to study the equilibrium features of this Lienard system, shows the existence of a saddle point in real time, with 
the implication that the stationary background solutions will be unstable, if the perturbation is extended into the nonlinear regime. 

So to summarise the import of this work, conservative momentum balance and continuity conditions, as appropriate for a stationary 
spherically symmetric flow, have been subjected to time-dependent radial perturbations. On including nonlinearity, an instability is seen 
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to develop in this otherwise simple hydrodynamic system. The entire mathematical treatment so described, and all its attendant physical 
conclusions, have been presented in what follows. 



2 THE MATHEMATICAL CONDITIONS OF SPHERICALLY SYMMETRIC ACCRETION 



The mathematical problem th a t was set up by Bondi ( 1952) himself and that is now taken up as a starting model in accretion-related 
texts Jchakrabartilll99(l Il996l : iFrank et alJbooaX involves two coupled fields, the local flow velocity, v, and the local density, p, of the 



compressible accreting fluid. These two coupled fields are governed by the continuity equation, 

dp 



+ 



1 d ( '\ n 

^irApvr )=0, 



dt r 2 dr 
and the inviscid Euler equation. 

1 dP 



dv dv 
dt dr 



+ --5J- + S'(r) = 0, 
p or 



(1) 



(2) 



tailored as they are, according to the requirements of spherical symmetry. In the latter equation, the local pressure, P, is expressed in terms of 
p, by invoking a general polytropic prescription, P — kp 1 , in which 7, the polytropic exponent, varies over the range (limite d by isothermal 
and a diabatic conditions), 1^7^ cp/cv, with cp and cv being the two coefficients of specific heat capacity of a gas jchandrasekhaj 
1939). The polytropic prescription is of a much more general scope than the simple conserved adiabatic case, and is suited well for the study 
of open systems like astrophysical flows. Now, making use of both P and p, it is also expedient to scale the flow velocity, v, in terms of a 
natural hydrodynamic scale of speed, c s , which is the local speed of sound. This speed can be noted from = dP/dp = jkp"' -1 . 

The flow is driven by the gravity of a central accretor, whose potential is $(r). In equation lO the driving force arising due to this 
potential is implied by its spatial derivative (represented by the prime). In the case of stellar accretion, the flow is driven by the Newtonian 
potential, <3?(r) = —GMfr. On the other hand, quite often in studies of accretion onto a non-rotating black hole, it becomes convenient 
to dispense with the rigour of general relativity, and instead make use of a pseudo-Newtonia n potential that mimics the general relativistic 
effects of Schwarzschild space-time geometry in a Newtonian construct of space and time 1 Paczvhski & Wiit J 198dl : Nowak & Wagoner 



199ll : lArtemova et al. 1996; Das & Sarkar 2001). The choice of a particular form of the pseudo-Newtonian potential, however, does not 



affect overmuch the general arguments regarding the stability of the flow. 

With the functions, P and §(r), specified, equations (T) and (0 can give a complete description of the hydrodynamic flow in terms 
of the two fields, v(r,t) and p(r,t). From these dynamic variables, the steady solutions of the flow are obtained by making explicit time- 
dependence disappear, i.e. dv/dt — dp/dt — 0. The resul ting differential equations, invo lving full spatial derivatives only, can then be easily 
integrated to get the stationary global solutions of the flow (Bondi 8l952l : Fran k et al . 2002). A remarkable feature of these stationary solutions 
is that they remain invarian t under the trans formation v — ► —v, i.e. the mathematical problem of inflows (v < 0) and outflows (v > 0) is 
identical in the steady state dChoudhurilll999h . This invariance has some adverse impl ications for critical flows in accretion pr ocesses. Critical 
solutions pass through saddle points in the stationary phase portrait of the flow dRav & Bhattacharieel2002ljRov & Ravl2007l) , but generating 
a stationary solution throu gh a saddle point will be imp ossible by any physical means, because it calls for an infinite precision in the required 
outer boundary condition 1 Ray & Bhattachariee1l2002 ). Nevertheless, criticality is not a matter of doubt in accretion processes 1 Bondilll952 ; 
Garlickll979l : lshapiro & Teukolskvll983nT The key to resolving this paradox lies in considering explicit time-dependence in the flow, because 
of which, as one may note from equations and l[2j, the invariance under the transformation, v — > —v, breaks down. Obviously then, a 
choice of inflows (v < 0) or outflows (v > 0) has to be made at the very beginning (at t — 0, as it were), and solutions generated thereafter 
will be free of all the difficulties associated with the presence of a saddle point in the stationary flow. 



On imposing various boundary conditions on the stationary integral solutions, multiple classes of flow result dFrank et alj|2002f) . Of 
these, the one that attracts attention in accretion studies obeys the boundary conditions, v — ► as r — ► 00 (the outer boundary condition) 
and v > c s for small values of r. It is quite obvious that this solution is transoni c in nature, with its bulk flow velocity overcoming the local 
speed of sound at a particular point in space, r c , the critical radius of the flow dchakrabartilll99ol : lFrank et al.ll2002 : iRav & Bhattacharjed 
2002). For a flow driven simply by the Newtonian potential, there is only one such critical radius. With the choice of a pseudo-Newtonian 
potential, multiple values of r c could result, but practically s peaking there would be only one physically relevant critical point, through which 
an integral solution co uld pass and a ttain the transonic state jMandal et alJ 2007). 

It was argued by Bondi 1 1952h that among all the feasible stationary solutions by which a fluid element may reach the accretor, after 
having starte d under highly subsonic conditions on very large length scales, the actual trajectory chosen will be the one that is transonic in 
nature — the Bondi ( 1952) solution. This line of thinking was based on the criteria that with no restrictive inner bou ndary conditio n, the 
accretion rate will be as high as possible and the corresponding energy configuration of the flow shall be the lowest one jGarlickll 1979b . The 
transonic solution conforms to these requirements, tak ing into co nsideration only the stationary conditions. Under the approximation of a 
"pressureless" motio n of a fluid in a gravitational field (Shu 1991), qualified support for transonicity also came later from a non-perturbative 
dynamic perspective 1 Rav & Bhattacharieell2002 : Roy_&Ra^ 2007 ). No definitive conclusion about transonicity, however, can be drawn on 
the basis of a perturbative linear stability analysis iGarlickl *1979h . 



Now, so far as generating the transonic flow is concerned, the non-perturbative dynamic evolution of global v(r, t) and p(r, t) profiles 
is very crucial indeed. Certainly, all the feasible stationary inflow solutions obey the outer boundary conditions that on large spatial scales. 
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and p(r) — > poo, where p^ is the constant "ambient" value of the density field very far away from the accretor {Frank et al 



2002). It is the way in which the two fields evolve close to the accretor that determines if the transonic state would be achieved or not. The 



dynamic process should be envisaged mathematically as one in which both the velocity and density fields, v(r, t) and p(r, t), are uniform 
initially for all values of r, in the absence of any driving force. Then with the introduction of a gravitational field (made effective at t = 0), 
the hydrodynamic fields, v and p, start evolving in time. If the temporal growth of v outpaces the temporal growth of p (to which c s is 
connected) at small values of r, then the final st ationary infall proces s will be transonic. Otherwise, the final stationary infall process will be 
globally subsonic, with v(r) — > as r — > |PettersonetalJl9801) . 

The non-perturbative evolution of the velocity and density fields in spherically symmetric accretion, however, requires working with a 
coupled set of nonlinear partial differential equations, as implied by equations (TJ and And where nonlinear equations are involved, one 
has to tread with caution, especially since no analytical solution of the dynamic problem exists in the case of spherically symmetric accretion. 



3 NONLINEARITY IN THE PERTURBATIVE ANALYSIS 



Equations (T) and ((2) are easy to integrate in their stationary limits, and the resulting velocity and density fiel ds, derived from thes e two 
equations, have only spatial profiles, v = vo(r) and p = pair). A standard practice in perturbative analysis l Petterson et al. 1980) is to 



apply small time-dependent, radial perturbations on the stationary profiles, Voir) and po{r), and then linearise the perturbed quantities. This, 
however, does not offer much insight into the time-dependent evolutionary aspects of the hydrodynamic flow. So the next logical step is 
to incorporate nonlinearity in the perturbative method. With the inclusion of nonlinearity in progressively higher orders, the perturbative 
analysis incrementally approaches the actual time-dependent evolution of the global solutions, after it has started with a given stationary 
profile at t — (to make physical sense, this initial profile has to be very much subsonic at all spatial points). 

The prescription for the perturbation is v(r, t) = Vo(r) + v'(r, t) and p(r, t) = po(r) + p'(r, t), in which the primed quantities indicate 
a perturbation about a st ationary background. I t is n ow necessary to define a new variable, f(r, t) = pvr 2 , following a similar mathematical 



procedure employed by IPetterson et al. ( 19800 and Theuns & David 1 1992 ). This variable emerges as a constant of the motion from the 



stationary limit of equation (Q}. This constant, /o, can be identified with the matter flow rate, within a geometrical factor dFrank et"al1l2002h . 
and in terms of vo and po, it is given as fo = povor 2 . On applying the perturbation scheme for v and p, the perturbation in /, without losing 
anything of nonlinearity, is derived as 

/' P' v' ft V 

4~ = — + — + . (3) 

Jo Po vo po vo 

The foregoing relation connects the perturbed quantities, v' , p' and /', to one another. To get a relation between only p' and /', one has to 
go back to equation (TJ, and apply the perturbation scheme on it. This will result in 

2l = _±df 4 

dt r 2 dr 

To obtain a similar relationship solely between v' and /', one needs to combine the conditions given in equations ((3) and to get 

dv' v fdf , df'\ ... 

u-5- • (5) 



dt f \dt dr 

In equations ®, © and O, all orders of nonlinearity have been maintained. Adhering to the same principle, applying the perturbation 
scheme in equation {2} and taking its second-order partial time derivative will yield 

d 2 v' d ( dv' c 2 B d P '\ „ 

H ( v——— -j- — ^— ] = 0. (6) 



dt 2 dr \ dt p dt 

In deriving this expression, all the terms involved in the stationary flow have vanished due to taking a partial time derivative. This is slightly 
different from the practice of extracting the stationary part of equation lO and making it disappear by setting its value as zero. Now making 
use of equations l[4}, O and the second partial time derivative of equation {5]l, a fully nonlinear equation of the perturbation is obtained from 
equation l|6}, in a symmetric form going as 

d f ^T) + 1 (h^-f) + #- (h"°£) + 1 ( h - d -f) = o, (7) 



dt \ dt J dt \ dr J dr \ dt J dr \ dr 
in which, 

2 

,tt V tr rt V rr V , 2 2\ , , 

h = j, h = h = h = - [v - Cs) . (8) 
Going by the symmetry of equation Q, it can be recast in a compact form as 

d^ (h^d v f) = 0, (9) 

with the Greek indices running from to 1, under the equivalence that stands for t and 1 stands for r. Equation l(9}, or equivalently, equa- 
tion I0, is a nonlinear equation containing arbitrary orders of nonlinearity in the perturbative expansion. All of the nonlinearity is carried in 
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the m etric elements, involving the exact fi eld variables, v, c s and /, as opposed to containing only their stationary background counter- 
parts Jvisserll998lJSchutzhold & Unrubl2002h . This is going into the realm of nonlinearity, because v and c s depend on /, while / is related 
to /' . If one were to have worked with a linearised equation only, then W v could be read simply from the symmetric matrix dRov & Ray 
2007h . 



1 

vo 



vo 



vl 



2 
C s0 



(10) 



in which c s o (r) is the stationary value of the local speed of sound. Now, in Lorentzian geometry the d' Alembertian for a scalar field in curved 
space is expressed in terms of the metric, as 

1 



Aip = 



(ii) 



with g^ v being the inverse of the matrix implied by g M „ i Visserl 19981 ; Barcelo et al.l2005 ). Comparing equations ((9) and i ll It with each other, 
one could look for an equivalence between and yj— g g^ v . What can easily be appreciated from this comparison is that equation (|9j gives 
an expression for /' that is of the type given by equation d 1 1 b . In the linear order, the metrical part of equation l|9}, as equation d 1 Ob shows it, 
may then be extracted, and its inverse will incorporate the notion of the sonic horizon of an acoustic black hole, when v% — c^o . This point of 
view has features that are similar to the metric of a wave equation obtained by setting the velocity of an ir rotational, inviscid and barotropi c 
fluid flow as the gradient of a scalar potential, and then by imposing a perturbation on this scalar potential I Vissejl998 ; Barcelo et alj2005 ). 
In contrast to this approach of exploiting the conservative nature of the flow to craft a scalar potential, the derivation of equation (|9j makes 
use of the continuity condition. The latter method is more robust because the continuity condition is based on matter conservation, which is 
a firmer conservation principle than that of energy conservation, on which the conventional scalar-potential approach is founded. 

Either way, all of this indicates that the physics of supersonic acoustic flows closely corresponds to many features of black hole physics. 
All infalling matter crosses the event horizon of a black hole maximally, i.e. at the greatest possible speed. By analogy the same thing may be 
said of matter crossing the sonic horizon in spherically symmetric inflows. Indeed, a long-standing conjecture about spherically symmetri c 
accretion on to a point sink is that the transonic solution crosses the sonic horizon at the greatest possible rate i Bondilll952 : Garlick|[l979h . 
That this fact can be appreciated for the accretion probl em through a p erturbative result is remarkable, because conventional wisdom would 
have it that perturbative techniques are inadequate here l Garlickfl979h . 

However, all of this is valid only as far as the linear ordering goes. When nonlinearity is to be accounted for, then instead of equation d 10b . 
it will be equations ^ which will define the elements, h? v , depending on the order of nonlinearity that one wishes to retain (in principle one 
could go up to any arbitrary order). The first serious consequence of including nonlinearity is to lose the argument in favour of the transonic 
condition (an inflow solution crossing the sonic horizon), because the descrip tion of as stated in equation d 10b . will not suffice any 
longer. This view is in perfect conformity with a numerical study conducted bv lMach & M alec (2008) for the case of spherically symmetric 
accretion, in which it was shown that if the perturbations were to become strong then the analogy between the "sonic horizon" and the event 
horizon of a black hole would not hold. Nevertheless, a most remarkable fact has emerged in consequence of including nonlinearity in the 
perturbative analysis. It is that regardless of the order of nonlinearity that one may desire to go up to, the symmetric form of the Lorentzian 
metric equation will remain unchanged, as shown very clearly by equation l(9}. For the laborat ory fluid problem of the hydraulic jump, a 
similar type of symmetry was shown to exist, going up to the second order of nonlinearity (Rav & Bhattachariee 2007b). 



4 STANDING WAVES ON STEADY GLOBAL INFLOWS 

All physically relevant inflow solutions obey the outer boundary condition, v(r) — > as r — ► oo. In addition, if the solution is globally 
subsonic, then the inner boundary condition is v(r) — > as r — > 0. From the point of view of a gravity-driven evolution of an inflow 
solution to a transonic state, the subsonic flows have great importance, because the initial state of an evolution, as well as the intermediate 
states in the march towards transonicity, should realistically be subsonic. So the stability of globally subsonic solutions must have a significant 
be aring on how a transon ic solution will develop eventually. Imposing an Eulerian perturbation on subsonic inflows, their stability was studied 
bv lPetterson et al. I Jl980h . and the amplitude of the perturbation in this case was seen to maintain a constant profile in time. In that respect 



one may say that the solutions do not exhibit any obvious instability. However, it is never prudent to extend this argument too far, especially 
when one considers nonlinearity in the perturbative effects, as it rightly ought to be done in a fluid flow problem. 

Now equation ((7) gives a nonlinear equation of the perturbation, accommodating nonlinearity up to any desired order. T his equation can 
be app lied to study the stability of stationary subsonic flows in a nonlinear regime. Following the mathematical procedure of lPetterson et alj 



ipp 

SO) 



( 1980), the perturbation is designed to behave like a standing wave about a globally subsonic stationary solution, obeying the boundary 



condition that the spatial part of the perturbation vanishes at two radial points in the spherical geometry — one at a great distance from the 
accretor (the outer boundary), and the other very close to it (the inner boundary). 

The mathematical treatment involving nonlinearity is to be confined to the second order only (the lowest order of nonlinearity). Even 
simplified so, the entire procedure will still carry much of the complications associated with a nonlinear problem. The restriction of not going 
beyond the second order of nonlinearity implies that W" in equations l[8j will contain primed quantities in their first power only. Taken 
together with equation Q, this will preserve all terms which are nonlinear in the second order. So, carrying out the necessary expansion of 
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v — vo+v' , p — po + p' and f = fo + /' in equations ([8} up to the first order only, and defining a new set of metric elements, g M " — foW , 
one obtains 

d» (q^duf) = 0, (12) 

in which p and v are to be readjust as in equation ^9). In the preceding expression, the elements, q^ v , carry all the three perturbed quantities, 
p' , v' and /'. The next process to perform is to substitute both p' and v' in terms of /', since equation d!2t is over /' only. To make this 
substitution possible, first one has to make use of equation ([3} to represent v' in terms of p' and /' in all g M ". While doing so, the product 
term of p' and v' in equation ([3} is to be ignored, because including it will raise equation J 1 2b to the third order of nonlinearity. Once v' has 
been eliminated in this manner, one has to write p' in terms of /'. This can be done by invoking equation l[4j, with the reasoning that if p' and 
/' are both separable functions of space and time, with the time part being oscillatory (all of which are standard mathematical prescriptions 
in perturbative analysis), then 



f =*(r)f, 
Po Jo 



(13) 



with a being a function of r only (which lends a crucial advantage in simplifying much of the calculations to fo llow). The exact functional 
form of a(r) will be determined by the way the spatial part of /' is set up. It was shown bv lPetterson et al. I dl98(j) that cr(r) would indeed be 
a real function, going as a(r) = vo (vo ± c s o)~ , when the spatial part of /' was cast as a power series in the WKB approximation. In any 
case, it stands to reason that when both p' and v' are real fluctuations, a should likewise be real. 

Following all of these algebraic details, the elements, q MI/ , in equation J 1 2b . can finally be expressed entirely in terms of /' as 

'./ I tr 2 f -. . t tr f \ rt 2 ( -t . /-rt f | rr (2 2 \ 3 <■ rr f 



vo(l + ef 



(14) 



in all of which, e has been introduced as a nonlinear "switch" parameter to keep track of all the nonlinear terms. When e = 0, only linearity 
remains. In fact, in this limit one converges to the familiar linear result implied by equation dlOt . In the opposite extreme, when e = 1, 
in addition to the linear effects, the lowest order of nonlinearity (the second order) becomes activated in equation d!2t . and the linearised 
st ationary conditions of a "sonic horizon" get disturbed due to the nonlinear e-dependent terms. This very feature has been tested numerically 
bv lMach & Maled d2008l) . Equations d!4t also contain the factors, , all of which are to be read as 



r = e 



2a, e 



3 + ( 7 -2) J' 

v o 



(15) 



Taking equations ( 1 1 2b . < l 1 4b and J 1 5 b together, a nonlinear equation of the perturbation is obtained, completed up to the second order, without 
the loss of any relevant term. 

To render equation d!2t . along with all g M " and £ M ", into a workable form, it will first have to be written explicitly, and then divided 
throughout by vq. While doing so, the symmetry afforded by £ tr = £ rt is also to be exploited. The desirable form of the equation of the 
perturbation should be such that its leading term would be a second-order partial time derivative of /', with unity as its coefficient. To arrive 
at this form, an intermediate step will involve a division by 1 + t^ u (f'/fo), which, binomially, is the equivalent of a multiplication by 
1 — e£ (/'//o)> w i m a truncation applied thereafter. This is dictated by the simple principle that to keep only the second-order nonlinear 
terms, it will suffice to retain just those terms which carry e in its first power. The result of this entire exercise is 



dt 2 dr 



Vo ■ 



dt 



i d 

+ —IT 

vo or 
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2vo dr 



2 ■ 
- c s0 
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dr 
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dr 1 S 
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dt 



v Q dC df' 2 
2 dr dt 



r 3 8f 
dr 



2£ u f' 



vo ■ 



dt 



Or 



I 2 
^0 (^0 



2 ' 
c s0 



df_ 

dr 



= 0, 



(16) 



in which, if one were to se t e = 0, then what would remain would be the linear solution discussed in detail by 



Pcttcrson ct al 



198G) 



and Theuns & David i 1992 ). To progress further, a solution of /' (r, t ) , separable in space and time, is to be applied. This will bear the form, 
f'(r,t) — R(r)<f>(t). Using this separable solution in equation d!6t , then multiplying the resulting expression throughout by vqR, and then 
performing some algebraic simplifications by partial integrations, will finally lead to 



4>v R 2 + <j>^- {v R) 2 + <f> \ A 
dr dr 



fo 



V0 I 2 



2 x dR" 

C S ) -7— 

' dr 



( 2 
V0 (Vo 



2 \ 
C s0 ) 



dR 

dr 



h 2 £ u voR 3 + < 



d , e rt 2 D 3\ . t rtVo dR 3 ttr, d / , , 

Tr (Z v R)+Z y^r-e R Tr (voR) 



2 I I 2 2\ dR d , tt 2 \ /-rr 3 r> ( dR\ d 



tttVo , 2 2 \ dR 3 

t. T (Vo - C s0 ) -T- 



d f fi rr «0 dR 3 
d7 I Y~dV 



= 0, 



(17) 



q \"U "SU / J 

S v ' dr 

in which the overdots indicate full derivatives in time. Quite evidently, equation d!7t is a second-order nonlinear differential equation in both 
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space and time. The way forward now is to integrate all spatial dependence out of equation J17t . and then study the nonlinear features of 
the time-dependent part. The integration over the spatial part will necessitate invoking two boundary conditions, one at a small value of r 
(close to the accretor), and the other when r — > oo (v ery far from the accre tor). At both of these boundary points, the perturbation will 
have a vanishing amplitude in time. It was reasoned by Petterson et al. i 1980h that globally subsonic inflow solutions offer conditions for 
the fulfilment of the two required boundary conditions, and simultaneously maintain a continuity of the background solution in the interim 
region. The boundary conditions will ensure that all the "surface" terms of the integrals in equation (17) will vanish (which explains the 
tedious mathematical exercise to extract several such "surface" terms). So after carrying out the required integration on equation J 17b . over 
the entire region trapped between the two specified boundaries, all that will remain is the purely time-dependent part, having the form, 



+ e(A4> + B4>) (j> + C(l> + eVcjf 



0, 



(18) 



in which the constants, A, B, C and T> are to be read as 



A = 
B = 



C = - 
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v R dr 



v R z dr 



vqR dr 
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3 dr 
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e'Ri- (voRf 
dr 
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2 \ dR d . tt 2 



dr, 



t-TT 3 D / dR 



dr. 



(19) 



respectively. The form in which equation d!8t has been abstracted is that of a general Lienard system l StrogatJl994 ; Jordan & Smith! 1999h . 
All the terms of equation J 1 St . which carry the parameter, e , have arisen in consequence of nonlinearity. When one sets e = 0, one readily 



Petterson et al 



1980). However, to go beyond linearity, and to appreciate the role of nonlinearity 



regains the linear results presented by 
the perturbation, one now has to understand the Lienard system that equation d 1 8 b has brought forth. 



5 EQUILIBRIUM IN THE LIENARD SYSTEM 

The mathematical form of a Lienard system is like a damped nonlinear oscillator equation, going as JstrogatJl99l : |jordan & Smith!l999h 



4> + eH((f>,<j))<j> + V( ( t>) = 0, (20) 

in which, H is a nonlinear damping coefficient (the retention of the parameter, e, alongside H, attests to the nonlinearity), and V is the 
"potential" of the system (with the prime on it indicating its derivative with respect to (/>). In the present study, 

U{ct>,4>) = A<t> + B4>, (21) 
and 

V(0)=C^+eD^, (22) 

with the constant coefficients, A, B, C and T> having to be read from equations d 1 9b . 

To investigate the properties of the equilibrium points resulting from equation OOt , it will be necessary to decompose this second-order 
dif ferential equation into a coupled first-order dynamical system. To that end, on introducing a new variable, ip, equation j20\ can be recast 
i Jordan & Smith! 19991) 



as 



(j> = V 

ip = -e(A4> + Bi>) V> - (C(f) + eV4> 2 ) . (23) 

Equilibrium conditions are established with <j> = %j) = 0. For the dynamical system implied by equations i23\ . this will immediately lead to 
two equilibrium points on the <f)-ip phase plane. Labelling the equilibrium points with a ★ superscript, one can easily see that (<f>* , ip*) = (0, 0) 
in one case, whereas in the other case, (<j)*,ip*) — (—C/(eD), 0). In effect, both the equilibrium points lie on the line, = 0, and correspond 
to the turning points of V ((/>). Higher orders of nonlinearity will simply have the effect of proliferating equilibrium points on the line, ip — 0. 
For the present case of second-order nonlinearity, one of the equilibrium points is located at the origin of the (jj-tp phase plane, while the 
location of the other will depend both on the sign and the magnitude of C /T>. 

Having identified the position of the two equilibriums points, the next task would be to understand their stability. To do so, both 
equilibrium points are to be subjected to small perturbations, following which a linear stability analysis will have to be carried out. The 
perturbation scheme on both <j> and ip is cj> = <f>* + 8cj> and ip = ip* + Sip. Applying this scheme on equation {23}, and then linearising in 8(p 



© 0000 RAS, MNRAS 000, 000-000 



8 Sen and Ray 



and Stjj, will lead to the coupled linear dynamical system, 



- (5<P) = H 



-V" {<j>*)5<j> - eU(</)* ,ip*)5ip, 



in which V" (</>*) = C + 2eD<f>* . Using solutions of the type, 
Jacobian matrix of the dynamical system follow as 



2 



,H 2 



V"( 



(24) 

exp(oji) and 5tp ~ exp(ajt), in equations 124t . the eigenvalues of the 



(25) 



with = %((/>*, V 1 *) having to be evaluated at the equilibrium points. 

Once the eigenvalues have been determined, it is now a simple task to classify the stability of an equilibrium point by putting its 
coordinates in equation (25}. The equilibrium point at the origin has the coordinates, (0, 0). Using these coordinates in equation d25t . the two 
roots of the eigenvalues are obtained as uj = ±iyC. If C > 0, then the eigenvalu es will be purely imag inary quantities, and consequently, 
the equilibrium point at the origin of the (f>-tp plane will be a centre-type point jjordan & Smithlll999l) . And indeed, when the stationary 



spherically symmetric inflow solution, about which the perturbat ion is constrained to behave like a standing wave, is globally subsonic, then 

„2 



C > 0, because in this situation, vq < c^o JPetterson et al. 1980h . Therefore, the centre-type equilibrium point at the origin of the phase 
pl ane indicates that the s tanding wave will be purely oscillatory in time, with no change in its amplitude. This very conclusion was drawn 
bv lPetterson et al. | jl98Ch in their linearised analysis of the standing wave, and it could be arrived at equally correctly by setting e = (the 



linear condition) in equation J25l >. 

The centre-type point at the origin of the phase plane has confirmed the results known already. It is the second equilibrium point that 
offers some novelties. This equilibrium point is entirely an outcome of taking nonlinearity to its lowest order (the second order) in the 
standing wave. The coordinates of this equilibrium point in the phase plane are (— C/(eD), 0), and using these coordinates in equation (25}, 
the eigenvalues become specified as 



AC 
2V 



MY 
2V ) 



+ C. 



(26) 



Noting as before, that C > 0, and that A, C and T> are all real quantities, the inescapable conclusion is that the e igenvalues, uj, are real quan- 
tities, with opposite signs. In other words, the second equilibrium point is a saddle point jjordan & Smitrll999l) , and as such its implications 
may be far-reaching when it comes to generating the transonic solution. 

To understand this, the first thing to note is that if the magnitude of the temporal part of the perturbation exceeds a certain critical 
value, i.e. if \cj>\ > \C/T>\, then the perturbation will undergo a divergence in one of its modes. In other words, the stationary subsonic global 
background solution will become unstable under the influence of the perturbation. This is how it must happen in the vicinity of a s addle 
point, and higher orders of nonlinearity (starting with the third order) will not smother this effect I StrogatJl994 ; Jordan & Smith! 1999h . The 
best that one may hope for is that the instability may grow in time till it reaches a satu ration level imposed by a higher order of n onlinearity, 
a feature that has a precedence in the laboratory fluid problem of the hydraulic jump dVoloviklEooi : iRav & Bhattacharie ell2007bl) . 

While all of this gives the perturbative perspective, the implications of the saddle point for the non-perturbative evolutionary dynamics 
are also noteworthy. It is evident that there can be no transonic solution without gravity driving the infall process. So from a dynamic point of 
view, gravity starts the evolution towards the transonic state from an initial (and arguably nearly uniform) subsonic state, far away from the 
critical conditions for transonicity. If, however, the subsonic states are to encounter a sad dle point in th e real-time dynamics, then that should 
hold adverse implications for reaching a stable and stationary transonic end, which is the Bondi ( 1952) solution. 

To ponder on a final point regarding the Lienard system, under linearised conditions, the perturbation on globally subsonic flows 
maintains a constant amplitude. Viewed in the phase portrait, this feature translates into close d phase trajectories around a centre- type point. 
Now, from dynamical systems theory, centre-type points are known to be "borderline" cases JStrogatj[l994l : |jordan & Smitlj|l9991) . In such 
situations, th e linearised treatment will show ap parently stable behaviour but an instability may emerge immediately on accounting for 
nonlinearity jStrogatJ 19941: 1 Jordan & Smitljfl 999). This is exactly what has happened in the perturbative study carried out here. 



6 CONCLUDING REMARKS 

Going by the form of the Lienard system derived in this work, it is easy to see that the number of equilibrium points will depend on the order 
of nonlinearity that one may wish to retain in the equation of the perturbation. In practice, however, the analytical task becomes formidable 
with the inclusion of every higher order of nonlinearity. Going up to the second order, an instability in real time appears undeniable, but then 
one must realise that this conclusion has been made regarding a purely inviscid and conservative flow. Now, actual fluid flows have viscosity 
as another important physical factor to influence their dynamics. In fact, fluid flows are usually affected both by nonlinearity and viscosity, 
occasionally as competing effects, and apropos of this point, it is to be noted that for a linearised pert urbation in spherically symmetric 
inflows, viscosity helps in decaying the amplitude of the standing waves on globally subsonic solutions i Ray 20031) . So the instability that 



© 0000 RAS. MNRAS 000, 000-000 



Implications of nonlinearity for spherically symmetric accretion 9 



has been seen to arise because of nonlinearity could very well be offset by accounting for viscosity in the flow. This is not to say though 
that viscosity will always act as a savio ur to preserve stability, because in one of the proposed models of axisymmetric accretion, viscosity 
has been known to destabilise the flow iBhattachariee & Ray 2007 : Bhattachariee et ai1l2009f) . The contrasting role of viscosity goes much 
beyond questions of stability. Looking at the respective geome tries in sph erically symmetric flows and axisymmetric flows, one notices that 
while viscosity te nds to inhibit the infall process in the former dRavl 20031) , it aids infall in the latter dshakura & Sunvaevlll973 ; Pringlell98 



Frank et al 



turbulence ( Meszaro 



2002h. Ap art from viscous dissipation , the stability of accretion p rocesses can, moreover, be af fected by radiative proc esses 
sfll975l ; lMeszaros & Silkll977l : lRav & Bhattacharieell2005n and magnetohydiodynamics (Balb us & Hawlevll 1998b. 



As a matter of regular practice, stability of fluids is also studied by constraining a perturbation to behave like a travelling wave dPetterson et al 



1980l : ICrosslll986l : Ray & Bhattacharje c 2007b). At time s, one encounters the surpr i sing situatio n of a fluid flow being stable under one type 



of perturbation, but unstable under the effect of another l Cross & Hohenber3ll993 ; Ray & Bhattacharje c 2007b). With nonlinearity lending 
an additional aspect, these effects merit a close examination in future studies. 



ACKNOWLEDGEMENTS 

This research has made use of NASA's Astrophysics Data System. The work here was carried out as a part of the National Initiative on 
Undergraduate Science (NIUS), conducted by Homi Bhabha Centre for Science Education, Tata Institute of Fundamental Research, Mumbai, 
India. The authors express their gratitude to J. K. Bhattacharjee, Tapas K. Das and S. Roy Chowdhury for helpful comments, and to Anwesh 
Mazumdar for support in various academic matters. 



REFERENCES 

Artemova I. V., Bjornsson G., Novikov I. D., 1996, ApJ, 461, 565 
Axford W. I., Newman R. C, 1967, ApJ, 147, 230 

Babichev E., Dokuchaev V., Eroshenko Yu., 2004, Phys. Rev. Lett., 93, 021 102 
Balazs N. L., 1972, MNRAS, 160, 79 

Balbus S. A., Hawley J. F, 1998, Reviews of Modern Physics, 70, 1 
Barcelo C, Liberati S., Visser M., 201 1, Living Rev. Relativity, 14, 3 
Begelman M. C, 1978, A&A, 70, 53 

Bhattacharjee J. K., Bhattacharya A., Das T. K., Ray A. K., 2009, MNRAS, 398, 841 
Bhattacharjee J. K., Ray A. K., 2007, ApJ, 668, 409 
Blumenthal G. R., Mathews W. G, 1976, ApJ, 203, 714 

Bonazzola S., Falgarone E., Heyvaerts J., Perault M., Puget J. L., 1987, A&A, 172, 293 

Bonazzola S., Perault M., Puget J. L., Heyvaerts J., Falgarone E., Panis J. F, 1992, Journal of Fluid Mechanics, 245, 1 
Bondi H., 1952, MNRAS, 112, 195 
Brinkmann W., 1980, A&A, 85, 146 

Chakrabarti S. K., 1990, Theory of Transonic Astrophysical Flows, World Scientific, Singapore 
Chakrabarti S. K., 1996, Physics Reports, 266, 229 

Chandrasekhar S., 1939, An Introduction to the Study of Stellar Structure, The University of Chicago Press, Chicago 

Choudhuri A. R., 1999, The Physics of Fluids and Plasmas: An Introduction for Astrophysicists, Cambridge University Press, Cambridge 

Cowie L. L., Ostriker J. P., Stark A. A., 1978, ApJ, 226, 1041 

Cross M. C, 1986, Phys. Rev. Lett., 57, 2935 

Cross M. C, Hohenberg P. C, 1993, Reviews of Modern Physics, 65, 851 

Das T. K., 1999, MNRAS, 308, 201 

Das T. K., 2000, MNRAS, 318, 294 

Das T. K., 2004, Classical and Quantum Gravity, 21, 5253 

Das T. K., Bilic N., Dasgupta S., 2007, JCAP, 06, 009 

Das T. K., Sarkar A., 2001, A&A, 374, 1150 

Foglizzo T, 2002, A&A, 392, 353 

Frank J., King A., Raine D., 2002, Accretion Power in Astrophysics, Cambridge University Press, Cambridge 

Gaite J., 2006, A&A, 449, 861 

Garlick A. R., 1979, A&A, 73, 171 

Holzer T. E., Axford W. I., 1970, ARA&A, 8, 31 

Jordan D. W., Smith P., 1999, Nonlinear Ordinary Differential Equations, Oxford University Press, Oxford 
Kazhdan Y. M., Murzina M., 1994, MNRAS, 270, 351 
Kovalenko I. G, Eremin M. A., 1998, MNRAS, 298, 861 

Landau L. D., Lifshitz E. M., 1987, Course of Theoretical Physics - Fluid Mechanics, Butterworth-Heinemann, Oxford 
© 0000 RAS, MNRAS 000, 000-000 



10 Sen and Ray 



Mach P., Malec E., 2008, Phys. Rev. D, 78, 124016 
Malec E., 1999, Phys. Rev. D, 60, 104043 
Mandal I., Ray A. K., Das T. K., 2007, MNRAS, 378, 1400 
Markovic D., 1995, MNRAS, 277, 11 

McCray R., 1979, Active Galactic Nuclei (ed. Hazard C. & Mitton S.), Cambridge University Press, Cambridge 

Meszaros P., 1975, A&A, 44, 59 

Meszaros P., Silk J., 1977, A&A, 55, 289 

Michel F. C, 1972, Astrophys. Space Sci., 15, 153 

Moncrief V., 1980, ApJ, 235, 1038 

Naskar T., Chakravarty N., Bhattacharjee J. K., Ray A. K., 2007, Phys. Rev. D, 76, 123002 

Novikov I. D., Thorne K. S., 1973, Black Holes: Les Houches (ed. DeWitt C), Gordon & Breach, New York 

Nowak A. M., Wagoner R. V., 1991, ApJ, 378, 656 

Paczynski B., WiitaP. J., 1980, A&A, 88, 23 

Park K., Ricotti M., 201 1, ApJ, 739, 2 

Parker E. N., 1958, ApJ, 123, 664 

Parker E. N., 1966, ApJ, 143, 32 

Petterson J. A., Silk J., Ostriker J. P., 1980, MNRAS, 191, 571 
Pringle J. E., 1981, ARA&A, 19, 137 
Ray A. K., 2003, MNRAS, 344, 1085 

Ray A. K., Bhattacharjee J. K., 2002, Phys. Rev. E, 66, 066303 
Ray A. K., Bhattacharjee J. K., 2005, ApJ, 627, 368 

Ray A. K., Bhattacharjee J. K., 2007a, Classical and Quantum Gravity, 24, 1479 

Ray A. K., Bhattacharjee J. K., 2007b, Physics Letters A, 371, 241 

Roy N, 2007, MNRAS, 378, L34 

Roy N, 201 1, MNRAS, 413, 2873 

Roy N, Ray A. K., 2007, MNRAS, 380, 733 

Ruffert M., 1994, ApJ, 427, 342 

Salpeter E. E., 1964, ApJ, 140, 796 

Schutzhold R., Unruh W., 2002, Phys. Rev. D, 66, 044019 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Shapiro S. L., Teukolsky, S. A., 1983, Black Holes, White Dwarfs and Neutron Stars, Wiley, New York 
Shu E K., 1991, The Physics of Astrophysics, Vol. II: Gas Dynamics, University Science Books, California 
Silich S., Tenorio-Tagle G, Hueyotl-Zahuantitla E, 2008, ApJ, 686, 172 
Singha S. B„ Bhattacharjee J. K., Ray A. K., 2005, Eur. Phys. J. B, 48, 417 
Stellingwerf R. E, Buff J., 1978, ApJ, 221, 661 

Strogatz S. H., 1994, Nonlinear Dynamics and Chaos, Addison- Wesley Publishing Company, Reading, MA 
Theuns T., David M., 1992, ApJ, 384, 587 

Titarchuk L., Mastichiadis A., Kylafis N. D., 1996, A&A, 120, 171 
Titarchuk L., Mastichiadis A., Kylafis N. D., 1997, ApJ, 487, 834 

Toropin Yu. M., Toropina O. D., Savelyev V. V., Romanova M. M., Chechetkin V. M., Lovelace R. V. E., 1999, ApJ, 517, 906 
Tsuribe T„ Umemura M., Fukue J., 1995, PASJ, 47, 73 
Visser M., 1998, Classical and Quantum Gravity, 15, 1767 
Vitello P., 1984, ApJ, 284, 394 
Volovik G. E., 2005, JETP Lett., 82, 624 

Volovik G. E., 2006, Journal of Low Temperature Physics, 145, 337 

Wong K.-W, Irwin J. A., Yukita M., Million E. T, Mathews W. G, Bregman J. N, 201 1, ApJL, 736, L23 
Zampieri L., Miller J. C, Turolla R., 1996, MNRAS, 281, 1 183 

This paper has been typeset from a TpX/ LTjX file prepared by the author. 



© 0000 RAS, MNRAS 000, 000-000 



